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Abstract 

We study a chain of non-linear, interacting spins driven by a static and a time-dependent mag- 
netic field. The aim is to identify the conditions for the locally and temporally controlled spin 
switching. Analytical and full numerical calculations show the possibility of stochastic control if 
the underlying semi-classical dynamics is chaotic. This is achievable by tuning the external field 
parameters according to the method described in this paper. We show analytically for a finite spin 
chain that Arnold diffusion is the underlying mechanism for the present stochastic control. Quan- 
tum mechanically we consider the regime where the classical dynamics is regular or chaotic. For the 
latter we utilize the random matrix theory. The efficiency and the stability of the non-equilibrium 
quantum spin-states are quantified by the time-dependence of the Bargmann angle related to the 
geometric phases of the states. 

PACS numbers: 
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I. INTRODUCTION 



Advances in nanoscale fabrication of magnetic materials down to a finite chain of in- 
divide nnagnetie atoms Q triggered a nnnrber of studies on the ground state magnetic 
properties of finite, interacting spin chains For accessing the non-equilibrium states, 
in a linear chain one conventionally rotates the spins by applying a static magnetic field 
H and a variable magnetic field h(t) along a direction perpendicular to if ul- The spins 
affected by the fields are then deflected by an angle = U\T which can be desirably varied 
by changing the duration r of the field h(i). Here oj\ is the amplitude of h(t) in frequency 
units u)\ = /io7 (7 is the gyromagnetic ratio). For this scheme to be viable uj\ has to be in 
resonance with the system precessional frequency. In this paper we consider the spin deflec- 
tion in the different situation of a nonlinear chain of interacting spins HQ 

in which case the 

precessional frequency is dynamical and changes with the oscillation amplitude jqj. Hence, 
a control strategy 9] as in the linear chain case entails the use of chirped fields. Here 
we inspect a different route to spin control by exploiting the stochastic nature of the spin 
dynamics when appropriate fields are employed. This we show in a first step analytically. 
The advantage is that no special frequency tuning is used and more importantly the spin 



may be quasi stable at the deflected (non-equilibrium) angles when the 



iield h[t) is off which 



might be of interest for quantum information applications Q, Q y fl fl Q 

Disadvantage is the limited control of the switching time. Full numerical simulations con- 
firm our analytical predictions: Tuning the external fields such that the underlying classical 
spin dynamics is chaotic, stochastic switching occurs and a long-time quasi stabilization, 
i.e. a dynamical freezing (DF) of the deflected states is possible. Small fields cause only 
small fluctuations around the equilibrium state. For very strong fields, effects of magnetic 
anisotropy and exchange become subsidiary and hence the dynamics turn regular and no 
deflection with subsequent freezing is possible. 

For a finite spin chain we uncover analytically that our stochastic control (SC) scheme is 
governed by Arnold diffusion and give analytical expression for the Arnold diffusion coeffi- 
cient that in turn determine the time scale for SC. 

To inspect the influence of the quantum nature of the spins on our (classical) predictions 
we considered both the regular and the chaotic classical regimes and evaluated the so-called 
Bargmann angle which is a measure of the quantum distance between states in the Hilbert 
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FIG. 1: A schematics of the interacting spin chain. The magnetic anisotropy field sets the z 
direciton. Two magnetic fields are applied: a static {Hq )Z ) one along the z axis and a time- 
dependent field (Ho jX (t)) along the x axis. 

space and can be used to signal DF j^, . Using random matrix theory we prove indeed that 
SC and DF are possible at the driving field values that follows from our classical analysis. 



II. EQUATION OF MOTION 

A. Liouville equation for the spin chain 

Similar to the case studied in 4], |5|] we consider a system that can be modeled by a chain 
of N interacting (with coupling constant J) spin variables localized at sites j and having a 
uniaxial anisotropy with the anisotropy constant (3. Possible sources of the anisotropy field 
are discussed in Refs. [l| . Here we mention that the inclusion of a finite anisotropy is essential 
for the existence of a finite-temperature long-range order in the (infinitely long) chain. 
Further important consequences of the magnetic anisotropy for the phenomena discussed in 
this work are detailed below. The direction of the anisotropy field defines the z axis. A 
static (H Qz ) and a time-dependent (H 0x (t)) magnetic fields are applied along the z and x 
axes, respectively (cf. Fig{I]). H 0x (t) consists of Nk periodic pulses, i.e. 

H 0x (t) = eT^5(t - kT), (1) 

k 



where T is the period and eT is the field strength. As demonstrated explicitly [18( (for the 
classical spin dynamics) the shape (CE]) of the field mimics well the action of a finite-width 
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pulse as long as the pulse duration is smaller than the field-free precessional period of the 
spins. The time-integral over the field amplitude of the finite- width pulse sets the variables 
e Q. 

From the Hamilton operator [4], Q 

AT-1 N N N 

H J SjzSj+l,z + H 0x ^ Bjx + H 0z ^2 § i z + @ ^ ( 2 ) 

j=l j=l j=l j=l 

we find the spin equation of motion (EOM) to be 

d t sj = hj x sj, hj = (H 0x , 0, J[s j+ i z + Sj-u] + H 0z + 2f3s jz ) . (3) 

For large spins Sj with hsj(sj + 1) = (sy) and [H, sf] = we shift variables as (cf. Figfj], 
s? = 1) 



Six = s i± cos ip^ s iy = s i± _ sin <p it s ij± = yl-s? r (4) 

EOM for two spins N = 2 in term of the canonical action variables Si z , S2 Z and their 

conjugate angles cpi, <p 2 reads 

dV dV 
Si z = - £ -^—> <Pi = Ui(s jz ) + e- — , Ui(s jz ) = Js jz + 2(3s iz + H 0z , i,je {1,2}; 

N N +oo N +oo 

V := ^2Vi(s iz ,ipi) =^2v 0i (s iz ,Lpi) T ^ 5(t- kT) = ^2 s i± cos Lp, r T 8(t-kT). 

i=l i=l k=—oo i=l k=—oo 

(5) 

For V ^ the variables of actions are adiabatic invariants and hence are slow with respect 
to the angles typical time scale |f| . The idea now is to identify the regime of classical chaotic 



dynamics which we wil 
the Liouville equation 
i.e. 

.df 



do below. In this regime one may adopt a kinetic approach based on 
3,0] for the two-particle distribution function f(t, s± z , tpi, s 2z , tp 2 ), 



1 ~qI = (L + sLx)f, 
L = -iui{s 2z )^- ioJ 2 {si z ) 



d(pi dip 
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, dV d dV d \ _.( dV d dV d \ 
ds u dtpt dtpxdsu) \ds 2z dcp 2 dcp 2 ds 2z J' 

Eq. (6) is of a key importance for this study. Below we use the random phase approximation 
and some mathematical techniques to derive from eg. (6) the Fokker-Planck equation which 
allows to explore some chaotic features of the spin dynamics. 



B. Fokker-Planck formulation and the onset of the chaotic regime 



Expressing f(si z , cpx, s 2z , ^2) as a Fourier series over ifi and ip 2 we find 



f(s lz , <fl,S 2z ,<f 2 ) = J7^^2 fn,rn{s lz ,S 2z )e mipl e 



vmip 2 



(27T) 

t t 

fn,m(siz,s 2z ) = f n ,m(si z ,s 2z ) ■ exp[- in J ui{t')dt'] ■ exp [ - im J u 2 (t')dt'], 



uj^t') = wi(s 2e (0), w 2 (i') =w 2 (a u (t'))- (7) 

Hence, solution of the Liouville equation is cast formally as (the symbol (n', m'\Li(ti)\n, m) 
means the average over fast oscillating variables) 

fn',m'(siz,S 2z ,t) = f n ',m'(si z , S 2z ,0) — 

f <(„'-«) /^(fjdf i(m'-m)J^(t')dt' ,. f ... , , 

2^ / dtie e ( nm |Li(ti)|n,m)/ nim (si^s 2z ,tiJ. (8) 



-75 

n,m 







If the interaction energy with the variable field is small with respect to the other terms in 
eq.([2]) we can expand / in terms of the field strength e and account for leading terms only. 
The zero-order component has the form 

/o,o( s lz) s 2z, t) = /o,0 { s lz, S 2z , 0) — 
t t 1 t 1 

f -in f u;i(t')dt' -im f U2(f)df 

-ie2_^ dt x e ° e (0,0\L 1 (t 1 )\n,m)f n:m (s lz ,s 2z ,0)+ (9) 



n,m 



t ti *2 t 2 

in f wi(t')dt' im f u) 2 (t')dt' 



/f in J LJi[t')dt' im J w 2 {t')dt' 

dh dt 2 e *i e (0,0|L 1 (f 1 )|n,m)(n, 7^1^(^)10,0) / Ol o(su,s 22 ,0). 

n 







Now we write as a Fourier series taking the relevant frequency Q = ¥jL into account, 
Li{t) = J2 P Li tP exp(ipQt), L±- p = L* p . Inserting into (9) we find 

/o,o( s lz, S 2z , t) = fofl(si z , S 2z , 0) — 

t H tj 

v-^ /" -in J uj 1 (t')dt' -im j u) 2 {t')dt' 

-ie^jdtye ° e ° (0, 0|Li(*i)|n, m)f n , m (s lz , s 2z , 0) + (10) 



+ (-^) 2 



in / uii(t')dt' im f u) 2 (t')dt' 



/P in J u>i(t')dt' im J 0J 2 {V)dV 

dh / dt 2 e «i e «i (0,0|Li,(ti)|n,m)<n,m|Li_p(* 2 )|0,0) 

n,m,p 
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Introducing the notations 

Mti,t 2 ) = J ^(tyt', Mti,t 2 ) = J u 2 (t')dt' (ii) 

h h 

we write 

fo,o{siz,S2z,t) = fo,o{siz,S 2 z,0) 
> 

£ J rft 1 e-^ 1 ^' )e-^( tl '°)(0,0|L 1 (t 1 )|n,m)/ n , m ( Slz , S22 ,0) + (12) 



-%£ 

n,m 



+{-ie 

n,m,p 



t ti 

2 



J dh J dt 2 e-^(* 1 '^e- im ^ ( * 1 '* a 5(0,0|Li iP |fi J m)(ri,m|Li > _p|0 J 0) 



o o 



■e ipn{tl - t2) f , {siz,s 2z ,0). 
Averaging over the initial phases, i.e. F(s\ z , s 2z , r) = (^fo,o(si z , s 2z ,t)^j , and using the 

*2 

random phase approximation i/)i(ti,t 2 ) = / toi(t')dt' ~ <fi(t 2 ) — fi(t 2 ), we end up with 

h 

<^exp [in^ 1>2 {t 2 ,ti)]j \ = exp ( - (ti - t 2 )/r c ) exp ( - ma^^i - r 2 )). (13) 

Here r c is the correlation time of random phase. Taking eq.(12) into account we deduce 
then for the averaged two-particle distribution function F(t) the dynamical equation (up to 
a second order in the field strength e) 

r\ 771 

-iee-% £ e -(-i+- 2 )*(o, |Li(t) |n, m)f n , m (su, s 2z , 0) - 



9* 



n,m 
t h 



m,n,p Q 

(0, 0\L hp \n, m) (m, n\L x _ p |0, 0)F(si z , s 2 ^, r). 

The long time behaviour (r ^> r c ) is retrieved by shifting to the new variables r = t\ — t 2 ,t\ = 
t\ in (14) and integrating over r which yields 

— = -e 2 V 2 ■ —r (0,0\Li p \n,m)(m,n\L l _ p \0,0)F(s lz ,s 2z ,t). (15) 

or — + ilnui + muo — pit) 

n,m,p T c v 

For a further progress explicit expressions for the matrix elements 
(0, 0|Li iP |n, m)(m, n\Li- p \0, 0) are need. Following the standard procedure outlined 
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FIG. 2: The time evolution of s\ z and s 2z for the following parameters of the system: J = 0.2, 2(5 
0, 1, H 0z = 0.2, e = 0.057, = 100. f = 0.1, s lz (0) = 0.8, ^i(0) = 0, s 2z (0) = -0.8, <^ 2 (0) 
0, T = and K = K' = 0.45 • 10~ 4 . 



in 21| we find after some lengthy steps the following Fokker-Planck equation for 

F(s lz , s 2z ,t) 



dF(s lz ,s 2z ,t) 



S I 7^(1- 



dF <9 



'I- si. 



dF 



D 



e 2 7t e 2 T 



dt \ds lz y ~ ~ lz> ds lz ' ds 2z ^ 2z ds 2z V 2tt 



(16) 



Making the ansatz F(si z , s 2z ) = F 1 (si z )F 2 (s 2z ) the average values of the spin projections Sj Z 
is determined from 



d 

—s 



d , _ 2 ^ dFj(s jz 
^-jz J j jz Q t u ■ a „ V 1 s jz) a „ 



(9s- 3 (9s 



2Ds jz ; j = 1,2; 



si, 2z = si i2z (0)e 



-2Dt 



(17) 



As discussed in 22( (for the case without anisotropy field), essential for the validity of this 
diffusion type dynamics is that the underlying classical dynamics is chaotic in which case 
the above derivations are justified. 

The stroboscopic evolution of the spin variables before (to — t) and after (t + t) applying 
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the field pulses at t — to is expressed as [22 1 



s iz,n+i — s iz , n + eTJl - s iz ^ n siny^, 

Vi,n+1 = ¥i,n + ( JSj z , n +l + W s iz ~ H 0z )T - eT _f^g_ _ COS f>i >n - (18) 

1 — S °z,n 

which also set the 



The stability of the trajectories is deduced from the Jacobian matrix 
condition for the chaotic regime as (< ... > t means time average) 

|Ai| > 1, K>0, (19) 

(2 + K)± v/(A' + 2)2-4 (2 - K) ± y/(K - 2)2 - 4 
A i,2 - 2 ' 3 - 4 ~~ 2 ' 

if = K ^/(f3/J) 2 + (1 - 2(/3/J) 2 ) Sl± s 2± < cos<^ cos^ 2 > t , K = eT 2 J. (20) 

Hence we can tune to the chaotic regime by varying the external fields parameters, the 
constant of anisotropy /3, and the coupling constant J between adjacent spins. For evaluating 
averages of the form < ... > t averages over time correlation functions of the random phases 
< cos tpi cos <£>2 > should be considered. For the correlation term we proceed as follows: 
When deriving the diffusion equation we assumed that correlation time of random phases 
are small with respect to the diffusion scale r c << ^ = 1/Q. Taking into account that in 
this time scale values s\^± are slow in time, after averaging of correlation functions over the 
time interval At e (0, we obtain K = K ^/(p/J) 2 + (1 - 2(0/ J) 2 )§ 5 s 1± s 2 ±t c ( J). 



C. Discussions and numerical results 



Having discussed the analytical structure of the spin dynamics we compare the analytical 
predictions with full numerical simulations of the problem. Note, our system is such 
that the conditions of chaotic regime eq.(19),(20) can be realized for arbitrary 
small perturbation e > 0, Kq > 0. However, the smallest values of e, and 
corresponding K = K' that allows for an observable effect has to be found 
numerically. At first, the external fields are tuned to K = K' = 0.45 • 10~ 4 > 0. In 
accord with the analytical results stochastic switching of the initial spins direction occur 
accompanied with a subsequent long-time stabilization (see Fig. (2)). If the fields are such 
that Kq < K' (i.e. Kq is very small) switching does not happen (see Fig.flSD), i-e. s-i^z is 
still an adiabatic invariants; external fields lead to small fluctuations around the equilibrium 
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ft 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



t x 10 2 ps 

FIG. 3: Same quantities as in Fig.© however, J = 0.2, 2/3 = 0.1, H 0z = 0.2, e = 0.04, Q = 
100. f = 0.1, s lz (0) = 0.9, c^i(0) = 0, s 2z {0) = -0.9, <^ 2 (0) =0, T = % and K = 
0.32 • 10~ 4 < K'. 

state. We note that in Fig. ([3]) the anisotropy field is finite but its effect is hardly observable 
becomes of it smallness ((3/ J) 2 = 0.06. The regular (but non-integrable) regime is reached 
by applying very strong fields (eT ^> J, H oz ^> Js z > 2(5s z ) (cf. eq. ([3])). In this case no 
stochastic switching occurs (cf. FigJH). The eigenfrequency of the system is given by the 
constant magnetic field Uj(Si Z ) = JSi Z + H 0z + 2f3sj Z ~ H 0z . Physically, effects related to 
the exchange interaction and to the anisotropy field become negligible and we end up with 
the familiar resonant switching scheme (this is true only during the time when the external 
fields are on. Effects of exchange and anisotropy govern the subsequent field-free dynamics. 
A scheme for a field-induced deflection and freezing has been proposed in Ref . fis| ) . 



D. Finite spin chain 

The EOM for a finite spin chain governing the dynamics of each particular spin follows 
from ([2]) as 

ds iz _ JWj (ipi, s iz ) 
dt difi 
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/ x 10 2 ps 



FIG. 4: Same quantities as in FigfiS however, J = 0, = 0, H 0z = 0.2, e = 0.04, 9, = 100. f = 
0.1, s u (0) = 0.9, pi(0)=0, s 22 (0) = -0.9, ^ 2 (0) = 0, T = and K = 0. 

difi _ , . <9V* s te ) , 

— — — (x>j ^Si_i z , Sj 2 , Sj+ixJ + £ , [Zl) 

at os iz 
Ui (si- lz ,s iz ,s i+lz ) = Jsi-u + Js i+lz + 2(3s iz + H Qz , i = l,..N. N + 1 = N. 

If the variable field has a spatial extent such that only two spins in the chain, labeled 
(k, k + 1), are affected then we find 

ds iz _ OVi {(fi,s iz ) 
- - (ft,* + W) * ^ . 

d(pi , . . , dVi(ifi,s iz ) 

—— — LOi {Si-i z , Si z , Si + i z ) -+- [O^k -T <>i,k+l) £ ~ 



ere asj z 
s iz(£) — constant if k ^ i ^ k + l. (22) 

These equations show that the z component of the spins subjected to the pulse have to 
be determined self-consistently. The dynamics of the oscillation frequency of the spins 
transverse components (pi, i = 1...N is determined by the effective magnetic field as 

V> i (t)=^" = 7sH! ,, (t)- (23) 
10 
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FIG. 5: We consider a chain of 10 interacting spins. We show here the dynamics of the z components 
s$ >z and S6,z of the spins labelled 5 and 6. Other parameters are J = 0.2, 2/3 = 0, 3, Hq z = 
0.2, e = 0.04, ft = 100. f = 0.1, 4 n+1 (0) = 0.9, sf n (0) = 0.9, c^(0) = 0, T 
N = 10. 



1,N. 



Here 



H eff (t) 



1 

7s 



#o + J(S.. 



i—l,z 



S i+ x !Z ) 4 [2(3- (5 itk + ^, fc+ i) ^ ( '^ cos(^) ) Si. 



S, 



i± 



This indicates that the spins subject to the pulses exchange energy with their nearest neigh- 
bors (whose z components are nevertheless constant). This process depends on the values 
of the z components and on the effective frequency (t)s*(t); a demonstration of this 
phenomena is shown in FigfSl 

A further tool for controlling the diffusion process is to apply a constant field along the x 
axis. The system dynamics is then chaotic, even without the periodic series of pulses {4, 0]. 
Therefore, the z component of the spin is not an adiabatic invariant and the mechanism 
of dynamical freezing (DF) discussed above does not work. Figj6] illustrates that if the 
amplitude of the magnetic field applied along x axis is strong enough H 0x > H 0z then the 
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FIG. 6: The time dependence of the z component of § rmth spin for the parameters J = 0.2, Hq z 



0.2, (3 = H 0x = l, s z 2n+1 (0) = 0.9, sf„(0) = -0.9, <^(0) = 0, i = l,N, N = 10. 

longitudinal component of the spin performs fast oscillations. In the other opposite case 
H Qx < H 0z the orientation of the spin can be deflected but DF again is not possible (cf. Fig 
[ID- 
After deflection of the spin to a desired angle, one can completely freeze its orientation. 
Key point is the fact that, in the absence of pulses the z component of the spin projection 
is an integral of motion. The system of equations in this case has the form 

^ = — ^( S j-l,z + S j+l,z) s jy ~ H 0z Sj y , 

ds ■ 

= J {sj-i, z + s j+1>z ) s jx + H 0z s jx , (24) 



ds 



dt 



JZ -0. 
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FIG. 7: The same as in Fig. [6] however the parameters are changed to J = 0.2, Hq z = 0.2, H { 



O.r 



0.01, sf n+1 (0) = 0.9, 4 n (0) = -0.9, ¥>i(0)=0, i = l,N, N = 10. 



For simplicity we assumed here that (3 = 0. After solving(24) one obtains 



(t ) COS (u j0 t + (fo) . 

s jy = s jy{ t o) sin i^jot + tfo), 



Sjz Sjz (^o) : 



(25) 



where Uj(sj-u{Q)', Sj+u(o)^ = jfsj-u(O) + s j+lz (o)^ + H 0z = u j0 , and s jz (t ) corre- 
sponds to the desired orientation of the spin, achieved after the action of the pulses in the 
time interval < t < t (cf FigJEJ). 
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/ ps 



FIG. 8: The z components of the spins labelled 5 and 6 for the parameters: J = 0.2, Hq z = 
0.2, H 0x = 0.01, s|n+i(0) =0.9, sf n (0) = -0.9, ^(0) = 0, i=TJf, N = 10. 

E. Arnold Diffusion 

Results of previous section evidence that even in the case of a long spin chain the ori- 
entation of spins can be still controlled. This follows from resonance overlapping and the 
existence of diffusion. However, the question of what kind of diffusion we have is still out- 
standing. If the dimension of the system is more than N > 2, the dynamics is much more 
involved and the emergence of new physical phenomenon is expected. We recall the key idea 
of KAM theory: The size of the destroyed torus is small and the domain of their location is 
surrounded by invariant torus. This situation changes if an invariant torus crosses the do- 
main of the destroyed torus location. This is possible if and only if N > 2. The phenomenon 
of universal diffusion along the net formed due to the inter-tours crossing was discovered by 

n 

Arnold |6J. Here we consider the mechanism of the formation of the stochastic net in the 
case of a spin chain. 

H = H {s' 1 ,...s* N )+eV(<p 1 , <p N ). (26) 

Note, the frequencies of the unperturbed motion on the N dimensional torus is a function 
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of the three actions 



^(sJLi, sf , s* +1 ) = Jsf_i + Js z i+1 + 2/3s* + i — l,N. 



(27) 



We collect the resonant tours defined by the condition: 



N 



w = °' 



(28) 



where n-,- are integer numbers. For each set of numbers there exists a multitude of 
solutions s° = (si , ...s^). Each solution determines resonant tours. For the formation of 
the Arnold diffusion the absences of degeneracy is essential 

d 2 H 



det 



ds*dsj 



^0, i,j = l,N. 



(29) 



In the case of our system due to the form of the matrix 

/ 2(3 J . \ 

d 2 H 



ds^dsj 



J 2(3 J . 
J 2(3 J . 
J 2(3 . 



(30) 



/ 



the condition of the absence of a degeneration (29) leads to the polynomial expression 



det 



d 2 H, 



ds^dsj 



G(J,(3,N)^0. 



(31) 



The explicit form of the expression fl3T|) also depends on the system's size (in addition to the 
dependence on the parameters (3, J). For large systems we obtain the following asymptotic 
expressions 



G(J,P,N) = J N if J>{3, G{J,(3,N) = 2 N (3 N if J < p. 



N nN 



(32) 



From this relation we can conclude that the universal diffusion is possible for any nonzero 
J, (3, and identify the numeric results obtained for the spin chain with the Arnold diffusion. 
For an analytical estimation we consider the minimal possible dimension. Therefore, in what 
follows, without loss of generality we shall restrict ourselves by the case N = 3. From eq.(28) 
we find 

n\(jj\ + n 2 LU 2 + "-3^3 = 0, (33) 
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where each frequency depends on the three action according to (27). In the frequency space 
(u)i,U)<2, 003) eq. (33) determines a family of surfaces. On the energy surface we have 

110(4,4,4) =E. (34) 

This equation is also the equation determining the surface (33). Therefore, the resonant 
tours have a common parts along the curves, defined as the solutions of the set of equations 
(33), (34). The time dependent perturbation leads to a widening of these curves and to the 
formation of the stochastic net. 

In order to provide topological interpretation of this phenomenon we will consider simplest 
case of three spins. In this case the explicit form of Eq (33) and Eq (34) reads: 

(J + 2/3) (oof + uo 2 2 + wf) - 2^i^ 2 - - 2oo 2 uo 3 = E, (35) 

and 

niUJi + n 2 uj 2 + n 3 uj 3 = 0. (36) 

Here E = (8/3 2 + 4/3 J — 4J 2 )H — 3Hq z (2/3 — J), is a re-scaled energy. From the equation 
(36) one can exclude frequency oo\ = — ^-oo 2 — and rewrite the equation (35) as a 
function of the two frequencies (oo 3 ,oo 2 ) . Clearly, the shape of the implicit plot for oo 3 (oo 2 ) 
depends on the values of parameters (711,712,713). Therefore, we expect that the 003(002), 
plotted for different inner resonances (711,712,713) should cross at some points. Now one 
can construct implicit plots expressing frequencies as a function of each others for different 
resonances (cf. Fig. 9). We see that, in some points trajectories cross each other. Due to 
topological reasons, such a nodal points are possible if and only if systems dimension is at 
least N = 3 or higher. Nodal points are crossing points between different resonances. If 
an external adiabatic perturbation is applied the dynamic near the nodal points becomes 
unpredictable and this leads to the Arnold diffusion [6]. 

In the general dimensional case (N » 1), the geometrical interpretation is less illustra- 
tive and much more complicated. Since we have to deal with N — l dimensional hyper-curves 
in the N dimensional hyper-space, the basic concept is however the same [20]). This conclu- 
sion manifests fundamental features of the multidimensional nonlinear dynamical systems. 
The diffusive motion of the system in the stochastic net is named as Arnold diffusion. 
Therefore, the diffusion equation (16) is still justified. However, the coefficient of diffusion 
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CO, 




FIG. 9: Topological structure of the inner resonances on the frequency plane (0^3,^2) plotted for 
different resonances: (n-2 = ni,n^ = n{), (77,2 = 2n\,nz = 3ni), (722 = 3ni,ri3 = ni), (ri2 = 
5ni, 77,3 = 7ni). Further values: J = 0.2, 2f3 = 0.3. 

for Arnold diffusion is defined by an expression other than eq. (16), namely Q|: 

D A = E 2 eH 0z e-^ a(N) . (37) 

Here E is the system's energy, a(N) is a dimensional dependent scaling constant with an 
upper limit determined by the Arnold inequality relation [6] 

"W < 6jV(Ar _ 1) 2 +3 iV + 14 - (38 » 
Obviously, for N » 1 , a(N) 1— > and the coefficient of the Arnold diffusion takes the more 
simple form 

D A = ^eE 2 H 0z . (39) 

Comparing eq. (39) with the diffusion coefficient obtained for the case of two spins, i.e. eq. 
(16) we find 
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This relation is important in that it delivers information on when the mechanism of 
stochastic switching and dynamical freezing are more efficient > 1 for a long spin chain, 
as compared to the case of a pair of spins. 

F. Role of anisotropy field 

Here we discuss the connection between the anisotropy field and Arnold diffusion. For the 
Arnold diffusion to occur the Jacobi matrix has to be none-degenerate. We note however, 
that the Jacobi matrix becomes degenerate in some cases if the anisotropy field is absent, 
as can be inferred from the structure of the Jacobi matrix. For example in the simplest case 
of three spins , 



det 



d 2 H 



dSfdS* 



-AJ 2 (3 + 8p' 3 ^ 0, if p^O. 



Evaluating the determinant for different number of spins we find that with the anisotropy 
field being applied it is always non-degenerated, while for particular N, it becomes degen- 
erated in absence of the anisotropy field. 

III. LOW TEMPERATURE LIMIT 

In this section we consider the dynamics in the continuous limit. By a proper choice of 
pulse parameters we were able to deflect the spin orientation diffusively to a desired angle. 
Switching off the pulses, the dynamics remain quasi frozen (is equivalent to the spin z com- 
ponent being an integral of motion). The question we pose here is that what happens if upon 
stochastic switching and freezing we apply a constant magnetic field along x axis. We recall 
that applying a constant field to the equilibrium (initial) state along the x axis invalidates 
the use of the KAM theory and the mechanism of SC does not work (cf.eq. (j5]l ). Before 
we deal with this problem in more details we set the limits of continuous approximation. 
Due to the constant magnetic field, applied along the x axis, the s z component is not an 
integral of motion any more. Therefore, excitations similar to spin waves propagate along 
the spin chain. These waves are not completely analogous to spin waves because the non- 
conservation of s z (t) is not related to flip-flop processes, but requires a transversal magnetic 
field. However, the wavelength of such excitations can be evaluated in a manner similar to 



the spin waves case 



241 ] . If the wavelength is larger than the distance between the spins 
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FIG. 10: The longitudinal s z , and the transversal s± = ys%. + Sy spin components as a function 
of time. J = 0.05 , H$ z = 1, Hq x = 0.1. Graph a) shows the analytical solutions given by eq.(52), 
whereas in graph b) the numerical integration of the system of equations (|47p is depicted. 

A >> a a continuous treatment is justified. Taking into account the expression for the wave 
frequency 

A\ T\.Q Ott 

(41) 



4|J|5 2tt 
lo = — - — • —— • a, 



h A 

One concludes that the validity of the continuous approximation depends on the temperature 



Jh 

T « —. 

Kb 



(42) 



Here fee is the Boltzmann constant. Thus, the continuous approximation corresponds to a 
low temperature approximation. For anti- ferromagnetic materials FeCl2, or C0CI2 we have 
J = 1.23 ■ 10 12 Hz, from (42) we infer for temperature regime of the continuous model 
T < 3K. 
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FIG. 11: The same as in Fig. 10 with the same meaning of the labels. The parameters of the system 
are however changed to J = 0.05 ,Hq z = 0.2, Hq x = 0.1. 

Returning back to the spin chain in a static magnetic field along x axis, the EOM read 



Considering that 



J\pj—l,z ^ 'jf+l,z) ^ 'jy Ho z Sjy, 

dsi 



dt 



J\Sj-l,z + Sj+i^JSja; + HQ z Sj x — Ho x Sj z , (43) 



J2 — Ff e 



Sjy ► Sy^X, t), (44) 
Sy 2 ► S 2 (x, t) 



ds z (x,t) a 2 d 2 s 



■r.t 



ds z {x,t) a 2 d 2 s Xyt 
s j+ i; Z = s x (x, t) + a dx + y-g^". (45) 
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from (43) we deduce that 

ds. r „ / d 2 s 



- W = -Js y \2s M + ^)-H 0m s yi 
ds T „ / <9 2 s 



ds z 
dt 



( d 2 s z \ 
Js x {2s z + -^-a 2 j-H 0x s z + H 0z s x , (46) 



H0 X Sy. 



In the low temperature regime we can neglect quadratic terms in (46) and obtain 

x = -Sy - yyz, 

y = 5x — z + jxz, (47) 
z = y. 

Here we introduced the following notations 
Hq z 2J 

o = -jj—, 7 = T7 — ' t h 0xh s x = x, s y = y, s z = z. 

Eq.(47) is derived in the absence of anisotropy field. However for the role of the 
anisotropy field we remark the following: The Zeeman field applied along the z 
axis is very strong H§ > J |>S'f_ 1 | + J \ Sf +1 \ + 2(3\Sf\. The eigenfrequency is constant 
and we have no effect of a dynamical shift u { (S^ , Sf +1 , Sf) = JSf^ +JS z i+1 +2pSf t + 
Hq w Hq. Inclusion of a finite anisotropy field leads to a rescaling of the small 
parameter 7 in in Eq. (47), i.e. 7 = — > Hence, we conclude that in this 

particular case the anisotropy field has no principal dynamical effect. Eq.(47) 
with rescaled parameter 7 = 2J + 213 is still valid in the presence of anisotropy field. 

When solving (47), we assume for the initial values the spin orientations achieved after 
action of pulses. In order to obtain analytical solutions we utilize the canonical perturbation 



theory (cf., e.g. [251] ). The parameter 7 is assumed to be small. First step is to rewrite (47) 
in a canonical form. This can be done using the following transformation. (For details, see 
Appendix) 

Xi = 25x — 2z, 

Vx = 2\y, (48) 
1 

Zl = -x + z, 



\ = VTT5 2 , 
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Equations (47) assume then the form 



#7 7<5 3 
xi = -Xyi + ^^i^i ^ -yWV^ 

yi = Xxi + 7 ^^ A 3~ - zizi - + (49) 



We seek a solution of (49) having the structure 

Xl = Cx\ + C 2 x? + C z xf + C 4 z< 4) 



Vl = Cy\ + C 2 y? + + CSt ] (50) 

Zl = Cz\ + C 2 4 2) + C 3 4 3) + C 4 z{ 4) . . . , 

and in addition we use the re-scaled time 

t = Ul + h 2 C 2 + h 3 C 3 + ...). (51) 
A 

With an accuracy up to third order in 7 we find the solution of (47) to be 

5 cos r 7(1 + 2S 2 ) sin 2 (r) _ 7 2 <5(9 + 4<5 2 ) cos(r) sin 2 (r) 
Sx ~ + 2A 4 " 16A7 ' 

75cos(r) sin(r) 7 2 (1 + 4<5 2 )(-5 sin(r) + 3 sin(3r)) 
s y - sm(r) — I , (52) 

cost 75sin 2 (r) 7 2 (— 1 + 45 2 ) cos(r) sin 2 (r) 
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A 2A 4 16A 

Fig. 10 demonstrate that these analytical solutions are in a good agreement with the exact 
numerical simulation of the system (I47p . Fig. 10 evidences that a constant magnetic field 
results in oscillations of spin's longitudinal component in a controlled manner. If the ampli- 
tude of the magnetic field is small, nonlinear effects become more important (cf Fig. 11). 



IV. QUANTUM MECHANICAL CONSIDERATION AND PROBLEM OF DF 

As stated above, the classical analysis is useful if the atoms in the chain have a large 
magnetic moment. In fact, for a finite chain of manganese (Mn) atoms [l| the classical 
approach proved to be adequate This situation changes however, for small spins where 
the dynamics becomes dominated by quantum effects. What is needed for our quantum 
consideration is the structure of the energy spectrum in the regime where the underlying 
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FIG. 12: Bargmann angle 6 sit) as a function of re-scaled time t = t/(Ho z + J), Hq z = J = 0, 2 
calculated from (59). 

classical dynamics is chaotic 4j. The point of interest here is that whether quantum effects 
invalidate the SC and in particular die dynamical freezing. To this end we use the concept 

28] in the way done in Refs.[s[y] to study DF. Let us consider 



of quantum geometry [26 



two quantum states and e^^i. The distance in Hilbert space between them can be 



characterized by the quantity 26 



27 



28) 



£>i(W 1; tf 2 ) = mmj^i - e itp ^ 2 \\. 

The minimal phase (p m is found by exterimizing — e l</, \I/2|| and noting that — \&2 
(\I>i|^2) 1/2 , which yields 

exp (i<p m ) 

Therefore, we writes for the distance Di 



(*1 


*2> 




#2>| 



(53) 



Di(*i,* 2 ) = V / 2-2|<$ 1 |^ 2 >|. 
For the same purpose as for Di, one may also use the Fubini-Study metric 

D 2 2 (^(t),vI/(0))=4(l-|<M/(t)|^(0)>| 2 ). 
In both approaches the key quantity is the so-called Bargmann angle 

coa6 B (t) = A(t) = |(tf(0)|tt(t))|. 
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27] 



(54) 



(55) 




FIG. 13: The time evolution of the Bargmann angle as calculated using eq.(62) and averaged 

over the random quantum phases. The parameters of the driving fields and the spin chains are the 
same as in Figj2l i.e. we are in the stochastic switching regime. 

In analogy to the classical deflection in quantum control problems, the Bargmann angle can 
be considered as the "quantum deflection". Essential for further progress is our (classical) 



?his calls for the use of 
29L hdl . To this end we 



finding that SC and DF occur in the classical chaotic regime, 
random matrix theory (RMT) to inspect the quantum dynamics 
write (|2J) as 

H(t)=H + V(t), (56) 
where Hq is time independent. Now we employ the established Floquet-operator method 



30| and the quantum map and infer for the Bargmann angle 



1 / N \ 

^(N+ cos [t{<p n -<p m )] J. (57) 



71 5^ 771 

Here (p n stand for the eigen phases of the Floquet operators. N is the Hilbert space dimen- 
sion. From this relation it is evident that starting at t — with two completely coherent 
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states, i.e. A(t = 0) = 1, de-coherence sets in for t ^ 0. To put the classical predictions 
of the previous section into a quantum perspective we note the following: In the classical 
regular regime, as identified above, the time-dependent perturbation V(t) acts adiabatically 
and does not alter the structure or the symmetry of the quantum spectrum. In this case we 
expect the Bargmann angle to be time periodic with typical quantum revivals. In contrast, 
in the classical chaotic regime, V(t) changes qualitative the quantum spectrum. 

Starting from the (classically) regular case we conclude that, if pair excitations are ne- 
glected, then the energy spectrum of the unperturbed part has the form 

E n = (n-1- N/2)H 0z + ^J(N - An + 3) + SE. (58) 

With this spectrum we infer for the Bargmann angle (eq.(57)) the expression 

1 / sin(iV+l/2)t-sin(t/2) V 1 / cosg/2) - cos(iV + l/2)t \ 2 
A[t) ~W\ 2sin(t/2) J + W\ 2sin(t/2) )' [ ) 

As evident from Fig. 11, if the underlying classical dynamics is regular, the time dependence 
of the Bargmann angle is periodical and the system is characterized by quantum revivals. 



To deal quantum mechanically with the classically chaotic regime we follow Ref. 30J and 



employ a Gaussian orthogonal ensemble 



30| 



P((p u ...(p N ) = Y[ (tp n - tp m ) exp ( - ^2<flj 

n>m ^ n=l ' 



N 

(60) 



i>i 

We note here that in general the distribution function (60) includes correlations between all 
N levels. If the number of correlated levels is n, then the n < N level-correlated distribution 
function reads 

P n {ipi, ...if n )= — - J P n (<fl, ■ ■ ■ <p N )d<p n +l ■ ■ ■ d(fi N . (61) 

The structure of the expression (57) suggests that the second-order correlated level distri- 
bution function ^(y^n; y 9 ™) is sufficient (each term in the sum contains two phases). Upon 
straightforward calculations we reduced P2(<p n ,<p m ) to 

mi ^Pn ) j 

with 

N 



fc=l 
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and 

M<P) = 1/2 ^fc(^)exp 



(2 n n\y/n) 

H n ({p) are Hermite polynomials. For (A 2 (t)^ p ^ we find (for N ^> 1 



r> rn m — 1 n — 1 



2 m-n m | 

m=l n=l 



/— TV m-1 . ^ 

(LT"(tV2)f - £ E E ^J tm """ '(C-'C 2 ^)) 2 (62) 

m=l n=l ^ 

where L™(t 2 /2) are Laguerre polynomials. The constant C we find from the condition 
(A 2 (0)) p = 1. Dynamical freezing (DF) means then a stabilization over time of the quantum 
distance jl], quantified by the Bargmann angle. To test for this situation we numerically 
solve for (62); a typical example is shown in Fig. ffTSl . These calculations are performed 
for parameters appropriate for the classically chaotic regime, e.g. those of FigfSJ The 
interpretation of Fig. ( 1131) is that SC drives the system diffusively to the target state 
which in this case is orthogonal to 1^(0)). The Bargmann angle is therefore deflected 
within a time to determined by the diffusion constant to the value Qb^d) ~ 7r/2. DF is 
then evidenced by a small variations of Ob^d) f° r t > tp. 



A. Conclusions 



For an exchange-coupled, non-linear spin chain and in the presence of a (uniaxial) mag- 
netic anisotropy and external driving fields, stochastic switching is possible if the field pa- 
rameters are chosen such that the underlying classical dynamics is chaotic. The switching 
mechanism is identified to be the Arnold-type diffusion. This we concluded analytically and 
substantiated by full numerical semi-classical and quantum calculations. We also inspected 
the possibility of dynamical freezing, i.e. stabilizing the target state beyond the switching 
time. 

Acknowledgment: The project is financially supported by the Georgian National Foun- 
dation (grants: GNSF/STO 7/4-197, GNSF/STO 7/4-179). The financial support by the 
Deutsche Forschungsgemeinschaft (DFG) through SFB 762 and though SPP 1285 is grate- 
fully acknowledged. 



26 



APPENDIX A 



The canonical Lyapunov system has the general structure 



x = -\y + X(x, y,z 1 ... z m ); 
y = \x + Y(x,y,zi . ..z m ); 



z 



- b s j + Z s (x, y, Z\ . . . z m ); 

3=1 

s — 1, 2, . . .rn; m — n — 2) 



We seek a reduction of this system of equations using the canonical transformation 

x = -Sy - yyz, 
y = 5x — z + 'jxz, 
z = y. 

We are interested in the linear part of eq. ( ]A3j) . i.e. 

x = -Sy, 
y = 5x — z, 

z = y, 

From the coefficients of the equation flA3j) we can construct the following matrix 



/ 



0-50 
5 0-1 
1 



\ 



We consider now the linear transformation 



fi = lux + 7i 2 y + 7i3^, 

6 = 721^ + 722?/ + 723^, 

6 = l3\x + 7 32 y + 733 z. 



In the new variables eq. (A3) is cast as 

dt 



Xi^i % = 1,2,3. 
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Taking eqs. (IA5p . (1A3|) into account we infer from eq. (1A6[) that 



(an - Xi)jn + a 2 i%2 + «3i7i2 = 0, 
aula + (a 2 2 - Ai)7; 2 + «327i3 = 0, 

Ol37il + «237i2 + (a 3 3 - \)li3 = 0. 

Equating the determinant to zero 



-Xi -5 
5 -Xi -1 
1 -\ 



0. 



we find 



X 1>2 = ±iX, A 3 = 0, A = Vl + 6 2 . 
According to (IA7I) this gives the following solutions for the matrix 

7n = - 7i2 = ~iX, 7i 3 = 1; 

721 = S, - 722 = -iX, 723 = 1; 

1 

731 = - 722 = 0, 733 = 1. 

So the canonical transformation (IA5B has the form 



The inverse transformation reads 



£i = — 5x — iXy + z\ 

£ 2 = Sx + zAy + z; 
1 

& = + 2. 



6 + 6 



2/ = + 6), 

Z ~ 2A 2 



Using eq. flAlOj) we obtain then for equations of motion in the variable £j 



6 = «Afi, 
6 = 0. 



28 



To reduce our system of equations to the canonical form, one more transformation is needed 



xi = -(6 + 6), 
yi =«(£l-6j), 

Zl = £3. 



(A14) 



Taking eqs. (lA10p , (lAllj) into account we obtain 



X\ = 25x — 2z; 
y\ = 2Ay; 



(A15) 



z 1 = -x + z. 





The inverse transformation reads 



5 ( xi 

i 

2) 
1 



x 



y = Yx yu 



x 2 



(-x 1 + 25 2 z 1 ). 



(A16) 



Using ( 1A15I) . ( 1A16I) we can finally rewrite the set of equations in the canonical form 



#7 7^ 3 

xi = -aj/i + 2^yi x i - -y^y^ 

. 7 S(S 2 - 1) 7 ^ 2 

yi = Axi + — x^ - —x 1 



2jP 2 
A 3 Zl ' 



(A17) 
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1 ^ 



4A3^ lXl " 2V yiZl - 
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